Thermodynamics, strange quark matter, and strange stars 
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Because of the mass density-dependence, an extra term should be added to the expression of 
pressure. However, it should not appear in that of energy according to both the general ensemble 
theory and basic thermodynamic principle. We give a detail derivation of the thermodynamics with 
density-dependent particle masses. With our recently determined quark mass scaling, we study 
strange quark matter in this new thermodynamic treatment, which still indicates a possible absolute 
stability as previously found. However, the density behavior of the sound velocity is opposite to 
the previous finding, but consistent with one of our recent publication. We have also studied the 
structure of strange stars using the obtained equation of state. 
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I. INTRODUCTION 

During the 10-plus years which have elapsed since Wit- 
ten's conjecture |l| that strange quark matter (SQM), 
rather than the normal nuclear matter, might be the 
true ground state of quantum chromodynamics (QCD), 
many theoretical and observational efforts have been 
made on the investigation of its properties and poten- 
tial astrophysical significance Q). Because of the well- 
known difficulty of QCD in the nonperturbative do- 
main, phenomcnological models reflecting the character- 
istic of strong interactions are widely used in the study of 
hadron, and many of them have been successfully applied 
to investigating the stability and properties of SQM. One 
of the most famous models is the MIT bag model with 
which Jaffe et. al [|| find that SQM is absolutely sta- 
ble around the normal nuclear density for a wide range 
of parameters. A vast number of further investigations 
HiH;@] are performed with fruitful results. A recent im- 
portant result is that young millisecond pulsars are most 
likely to be strange stars rather than neutron stars ||. 
Another alternative model is the mass-density-dependent 
model with which Chakrabarty et al. obtained signifi- 
cantly different results ||f[0]]- However, Benvenuto and 
Lugones [jll| pointed out that it is caused by the wrong 
thermodynamic treatment. They added an extra term 
to the expression of both pressure and energy, and got 
similar results as that in the bag model. A recent inves- 
tigation indicates a link of SQM to the study of quark 
condensates |fl2| while a more recent work has carefully 
studied the relation between the charge and critical den- 
sity of SQM ||. 

Latterly, we have demonstrated that th e p revious 
treatments have unreasonable vacuum limits p4| . In ad- 
dition to this problem, there exist another serious prob- 
lem, i.e., the zero pressure does not appear in the lowest 
energy state. In fact, there are two important problems 
in the quark mass-density-dependent model. One is how 
to determine the quark mass scaling. The other is how 
to deal with the thermodynamics with density-dependent 



particle masses self consistently. We have mainly concen- 
trated on the first problem in Ref. E| . The present paper 
will concentrate more on the second problem. We find 
that the extra term provided in Ref. |n| should indeed 
be appended to the expression of pressure. However, it 
should not appear in that of energy according to both the 
general ensemble theory and basic thermodynamic prin- 
ciple. After our modification, the zero pressure point 
appears exactly at the lowest energy state, and thus the 
thermodynamics with density-dependent particle masses 
becomes self-consistent, which leads to completely differ- 
ent density behaviour of the sound velocity in SQM and 
different structure of strange stars. 

We organize this paper as such. In the subsequent sec- 
tion, we give detail arguments on why the additional term 
in the pressure should not appear in the energy. The ther- 
modynamic expression needed later are all derived care- 
fully in this section. Then in Sec. HI, we apply the new 



thermodynamic formulas to investigating strange quark 
matter, and find that the density behaviour of the sound 
velocity is opposite to the previous treatment but consis- 
tent with our recent publication. On application of our 
obtained equation of state, we integrate the equations of 
steller structure for strange stars in Sec. IV, and find that 
strange stars are dimensionally smaller and less massive 
than the previous calculation if SQM is absolutely stable. 
Sec.0 is a short summary. 



II. THERMODYNAMICS WITH 
DENSITY-DEPENDENT PARTICLE MASSES 



Let us explore directly from the general ensemble the- 
ory what the expression of pressure and energy should 
look like if the particle masses are dependent on density. 

We express the density matrix as 
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where S is the partition function, (3 is the reverse tem- 
perature, Ni are the particle numbers, and /Xj are the cor- 
responding chemical potentials. The microscopic energy 
ENi.a is a function of the system volume V, the parti- 
cle masses to*, the particle numbers Ni, and the other 
quantum numbers a. The pressure of the system is 

p=- V f d-Ejy^cA c -ff(g Jff< , ot -y; tW jy i ) 
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is the thermodynamic potential density which is generally 
a function of the temperature T, the chemical potentials 
/ij, and the particle masses m^. If the particle masses 
have nothing to do with the baryon number density n b = 
N/(3V) (N is the total particle number), we simply get 



p = -n. 



(4) 



If the masses depend on density or volume, one should 
have 



p= -n + n b —. 



(5) 



This is just the right thing the authors have done in Ref. 
lOl. The authors derive it as such: 
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For canonical ensemble, the particle numbers Ni keep 
fixed. This derivation is thus obvious. However, it is not 
so obvious for grand canonical ensemble because the par- 
ticle number is not necessarily constant when the tem- 
perature T and chemical potentials fj,i unchanged. We 
will give a more convincing derivation a little later. 

The additional term is of crucial importance for pres- 
sure balance. Not as done in Ref. (TTJ , however, the extra 
term does not appear in the expression of energy. Now, 
let's calculate the statistic average for the energy: 
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is the average number for particle type i. Therefore, the 
energy density of the system is 



E = — = + 2^ 
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where rij is the number density of particle type i: 
N dfi 



V 
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It is clear from Eq. (|ll]) that only when Eq. ([|) holds 
can one get the Eq. (8) in Ref. Q. Therefore, we should 
not, as done in Ref. Q|, use that expression to calculate 
the energy density. Instead, we will calculate E directly 
from Eq. ( |Ti"|) in this paper. 

For more evident arguments, let's see the following 
derivation starting from the basic derivative relation for 
an open system: 



d(VE) = Td(VS) - PdV + ^ mdNi 



(13) 



where S is the entropy density of the system. 

Choosing T, V, and {N} as the independent macro- 
scopic state variables, the combined statement of the first 
and second laws of thermodynamics, Eq. (|l3|), can be ex- 
pressed as 



d(VA) = -VSdT - PdV + Vi dN i, 



(14) 



where A = E — TS is the Helmholtz free energy density 
by which we have 
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This is a general expression for pressure. In obtaining 
the third equality, we have used the chain relation 



-V W f ( {ni 



Ni 
V 



= y2 n 3l — f({ n i})-> (18) resulting integration: 

{Ni} j 3 



At zero temperature, the corresponding thermody- 
namic potential density can be obtained from Eq. ( |20) ) 
by first taking the limit T — > and then carrying out the 



where / is an arbitrary function. 

According to the basic relation between thermodynam- 
ics and statistics, we have 



A = Q + ^2 HiUi 



(19) 



where n is the thermodynamic potential density. For a 
free Fermi system, it is 
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We thus have from Eqs. (gj), (y), and Q: 

«i = ^(M?-m?) 3/2 , 
S = mimF(xi), 
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where is the degeneracy factor which is 6 for quarks where 
and 2 for electrons. In order to include the interaction be- 
tween particles, we regard the particle masses as density- 
dependent, namely 
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mi = mi | Tib = ^2 n j/3 ■ (22) is the ratio of the Fermi momentum to the mass that re- 

j J lated to particle type i. With the hyperbolic sine function 

- . sh-^x) = ln(x + Vx 2 + 1), the functions Ffa), Gfa), 

Because we have chosen T, 7, and {AT,} as independent and f ,. ^ defined &g 

state variables, the chemical potential pn should also be 

regarded as a function of T and {n^}, namely 



Hi = Hi (T, {n k }) ■ 



(23) 



So, the total derivative of Sl(T, {^k}, {«t/c}) with respect 
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One can see, from Eqs. (33) and (p4|), that an addi- 
tional term appears in the pressure expression, but not 
in the energy expression. We can specially confirm this 
result further as such. 



From Eq. (fL3[), one has an alternative general expres- 
Here we have used the Eq. (|12|) and the fact that gion for pressu 7e- 



drib/drij = 1/3. 

Substituting Eq. (Il9|) into Eq. (|17|) gives 
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According to the Pauli principle and the rclativistic 
energy-momentum relation = y/p 2 + m|, the energy 
(28)density of the system at zero temperature should be 
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,„ n *which, when the integration is carried out, is just the 
^ 'same as Eq. (|33|). 
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Because the entropy is also zero (or constant) at zero 
temperature, we can substitute Eq. (|4l| ) into Eq. (40), 
and accordingly get 
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drih drrii 
which leads to Eq. (|3) exactly. 



III. PROPERTIES OF STRANGE QUARK 
MATTER IN THE NEW THERMODYNAMIC 
TREATMENT 



Having derived in detail the thermodynamics with 
variable particle masses in the previous section, we now 
apply it to the investigation of strange quark matter. 

As usually done in the previous literature 
|,|Jl^Jll|,|l|,|l|,0, We assume the SQM to be a Fermi 
gas mixture of u, d, s quarks and electrons with chemical 
equilibrium maintained by the weak interactions: 



d,s <-> u - 



d,. 



Because of these reactions, the chemical potentials 
/ij (i = u, d, s, e) should satisfy 



Mu + Me 



(45) 
(46) 



For the bulk SQM in weak equilibrium, the previous 
investigations got a slightly positive charge ||. Our re- 
cent study jl3| demonstrates that negative charges could 
lower the critical density. However, too much negative 
charge can make it impossible to maintain flavor equi- 
librium. Therefore, the charge of SQM is not allowed to 
shift too far away from zero at both positive and negative 
directions. Therefore, one also has another two equations 
for a given baryon number density rib'- 



3 (yiv 



= n b , 

Tie = 0. 
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The first is the definition of baryon number density; the 
second is from the charge neutrality requirement, rii (i = 
u, d, s, e) is related to fii and m, by Eq. (]32]). 

Because the results from lattice calculations show that 
quark matter does not become asymptotically free soon 
after the phase transition (instead, it approaches the free 
gas equation of state very slowly), one should consider 
the strong interaction between quarks in a proper way. 



We do this by including the interaction effect within the 
variable quark masses. Because of the characteristic of 
the quark confinement and asymptotic freedom of the 
strong interaction, one can write down the simplest and 
most symmetric parametrization for the quark masses 
u,d,s) @: 



m q (q 
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where m q o is the corresponding quark current mass, z is 
a fixed exponential. Previously, z is regarded as 1. Our 
recent study jl4| indicates that it is more reasonable to 
take z = 1/3. The parameter D is usually determined 
by stability arguments, i.e., at zero pressure (P = 0), 
the energy per baryon, E/rib, is great than 930 MeV 
for two flavor quark matter in order not to contradict 
standard nuclear physics, but less than 930 MeV for three 
flavor symmetric quark matter so that SQM can have 
the possibility of absolute stability. Obviously, the rang 
of D determined by this method depends on different 
thermodynamic treatments. Within the thermodynamics 
derived in the preceding section, D is in the range (155 — 
171 MeV) 2 when taking z = 1/3. 

Because the light quark current masses are very small, 
their value uncertainties are not important. So we take 
the fixed central values m u o = 5 MeV and rrido = 10 
MeV in our calculation. The electron mass is very small 
(0.511 MeV). As for s quarks, we take 80 and 90 MeV, 
corresponding respectively to D 1 / 2 = 156 and 160 MeV. 

For a given rib, we solve for [ii (i = u, d, s, e) from Eqs. 
d45| ) — (^8|), and calculate the energy density and pressure 
of SQM respectively from Eqs. (33) and (|3~i| ) with the 
quark messes replace by Eq. (f49|). 

Firstly, we draw the configuration of the SQM for the 
parameter set m s0 — 80 MeV and D 1 / 2 = 156 MeV in 
Fig. [I]. At high densities, all of the u, d, and s quarks 
tend to become a triplicate. When the density becomes 
lower, d fraction will increases while s fraction decreases, 
and will become zero at a definite density which is called 
critical density in Ref. (ll| because SQM does not exist 
below that density. The u fraction is nearly unchanged. 
It in fact increases very slowly. To keep charge neutrality, 
the electron fraction will also increase. However, because 
of its very small mass, the electron fraction is so little that 
we multiply by one thousand to draw it in the figure. 

In Fig. we show the density dependence of the energy 
per baryon, VE/N — E/rib, vs baryon number density 
rib for the parameter set I: m s o = 80 MeV, D 1 / 2 = 156 
MeV, and II: m s0 = 90 MeV, D 1 ' 2 = 160 MeV. For 
the first parameter set, SQM is absolutely stable while 
for the second set it is nearly meta-stable. The point 
marked with a circle 'o' is the zero pressure point where 
the system pressure becomes zero. It can be clearly seen 
that the zero pressure points are exactly located at the 
lowest energy state. In fact, this is a basic requirement of 
thermodynamics because one can obtain from Eq. (|39| ) 
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However, this is not the case for most of the previous 
thermodynamic treatments of strange quark matter in 
the mass-density-dependent model ]ll|j9|, ^0| , p^ | , which is 
their another serious flaw in addition to the unreasonable 
vacuum limits mentioned before. 

In Fig. ||, we give the relation between the pressure P 
and energy density E, i.e., the equation of state. It ap- 
proaches the free gas equation of state at high densities. 
However, its shape is a little sunken at lower densities, 
contrary to previous calculation which is protuberant. 
This will leads to completely different lower density be- 
haviour of the sound velocity in strange quark matter. 

The velocity of sound is plotted in the lower part of 
Fig. [|. The upper part is calculated by the same method 
in Ref. JO] with parameter set B there. Simultaneously 
given with a full horizontal line is the ultra-relativistic 
case (1/V3) for purpose of comparison. Obviously, they 
become nearly identical at high densities while the lower 
density behaviour is opposite. The sound velocity in the 
previous treatment is higher than the ultra-relativistic 
case and will eventually exceed the speed of light at lower 
densities, which is unreasonable from the point of view 
of the theory of relativity. 



IV. STRANGE STARS 



It has long been proposed that the currently called 
neutron stars might be composed of strange quark mat- 
ter and thus be in fact strange stars. A recent investiga- 
tion shows that young millisecond pulsars are most likely 
to be strange stars rather than neutron stars Previ- 
ous authors have investigated the properties of strange 
stars by applying their obtained equation of state with 
interesting results pO|JTl] ] . We have now modified the 
thermodynamic treatment and updated the quark mass 
scaling. Therefore, it is meaningful to study the structure 
of strange stars in the new context from the astrophysical 
view point. 

As generally done, we assume the strange star to be a 
spherically symmetric object. Its stability is determined 
by the general relativistic equation of hydrostatic equi- 
librium known as Tolman-Oppenheimer-Volkov equation 
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dp 

dr 



GmE (1 + Anr 3 P/m)(l + P/E) 



r 2 1 — 2Gm/r 

with the subsidiary condition 

dm/dr = Att^E, 



(51) 
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where G = 6.707 * 10 45 MeV 2 is the gravitational con- 
stant, r is the distance from the core of the star, E = E(r) 



is the energy density or mass density, P — P(r) is the 
pressure, and m — m(r) is the mass within the radius r. 

For an initial baryon number density no (accordingly 
Pq and Eq), we can numerically solve Eqs. ( J5l| ) and (|52| ) 
with the aid of the equation of state, and obtain the 
corresponding P = P(r,no) and m = m(r, no), and con- 
sequently n = n(r, no), the baryon number density at the 
radius r for the central density no- The radius R of the 
strange star is determined by the condition 



namely, 



P(R,n ) = 0, 



R = R(n ). 



Accordingly, the mass of the strange star is 
M = m(R(n ),n ) = M{n ). 



(53) 



(54) 



(55) 



To make strange stars stable, we must require 
dM/driQ > 0. For the above obtained equation of state, 
M firstly increases with no up to a definite value M max 
corresponding to the highest acceptable central density 
nomax- After that, M decreases with no, and the star 
becomes unstable. 

For the parameter set I, i.e. ni s o = 80 MeV and 
D = (156 MeV) 2 , we give the density profiles n(r, no) in 
Fig. |5] as an example. The upmost line is for the largest 
acceptable central density no max (« 1.35 fm~ 3 ). The 
lowest horizontal line corresponds to the surface density 
n s (w 0.25 fm ~ 3 ) of strange stars which is independent 
of the central density, but a functional of the equation of 
state. Each line will intersect with it. The cross points 
correspond to the radius R of the star. The maximum 
radius of the star appears in no ~ 0.65 fm~ 3 . 

In Fig. ^, we show the mass-radius relation of strange 
stars with a solid line. The point marked with a full 
dot '•' represents the largest acceptable mass M max (« 
1.58 times the solar mass). For comparison, we have 
also plotted the result from the bag model calculation 
with the bag constant B X I A = 144 MeV, and that in 
Re. |ll[] with the parameter set B there. We can see 
that the strange stars in our case is dimensionally smaller 
and less massive than the previous calculation if SQM is 
absolutely stable. Naturally, this observation depends on 
the parameters employed. If we choose a bigger m s o and 
larger D, the case would be different. However, SQM 
would have no possibility of absolute stability in that 
case. 



V. SUMMARY 



We have derived the thermodynamics with density- 
dependent particle masses self-consistently, which over- 
come the serious flaws of the previous treatment of SQM 
in the quark mass-density-dependent model. We find 
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that an additional term should be appended to the ex- 
pression of pressure, but it should not appear in that of 
energy. When applying the new formulas to the inves- 
tigation of SQM, we find that the density behaviour of 
the sound velocity is opposite to the previous calculation, 
but consistent with our recent publication, which leads 
to different structure of strange stars. 
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FIG. 1. The configuration of SQM varies with density. At 
high densities, all of the u, d, and s quarks tends to become 
a triplicate. When the density becomes lower, d fraction will 
increases while the s fraction decreases. The u fraction is 
nearly unchanged. It only increases very slowly. 

FIG. 2. Energy per baryon E/rib vs baryon number den- 
sity rib for differnt parameter groups. The zero pressure 
points marked with 'o' are located at the lowest energy 
state, which is not the case for most of the previous ther- 
modynamic treatments of strange quark matter in the quark 
mass-density-dependent model. 

FIG. 3. Equation of state of strange quark matter (pressure 
P vs energy density E) . it approaches the free gas equation of 
state at high densities. However, it is a little sunken at lower 
densities, contrary to previous calculation. 
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FIG. 4. Velocity of sound in strange quark matter. The 
solid horizontal line is the ultra-relativistic case, the lower 
half part is the results of our calculation while the upper part 
is calculated by the same method in Re. jll| for parameter set 
B there. Their lower density behaviour is obviously opposite. 

FIG. 5. Density profiles for parameter set m s o = 80 MeV 
and D = (156 MeV) 2 . The upmost line is for the largest 
acceptable central density no mal . The lowest horizontal line 
corresponds to the surface density of strange stars. The cross 
points of each line and the lowest horizontal line correspond 
to the radius R of the star. 

FIG. 6. Mass-Radius relation for strange stars. The verti- 
cal axis is the star mass in unit of the solar mass while the 
horizontal axis is the star radius in unit of kilometer. The 
solid line is obtained by the method in this paper. The dot- 
ted line is from the bag model. The dashed line is calculated 
with the same method in Ref. Jll| for the parameter set B 
there. The points marked with a full dot '•' represent the 
maximum acceptable masses. 
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